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Abstract 


A  common  problem  in  the  field  of  enviromental  fluid 
mechanics  is  the  disposal  of  an  effluent  into  a  moving 
ambient  fluid.  The  prediction  of  concentration  contours  is 
important  in  these  cases  to  be  able  to  evaluate  the  impact 
of  the  effluent.  In  cases  where  the  effluent  discharge 
velocity  is  large  compared  to  the  ambient  velocity  the 
problem  may  be  modeled  as  a  turbulent  jet.  The  purpose  of 
this  study  is  to  investigate  the  effect  of  a  weak  crossflow 
on  such  a  jet . 

The  study  was  done  by  first  performing  an  experimental 
investigation.  The  situation  modeled  in  the  experiments  was 
that  of  a  jet  discharged  on  the  bed  of  a  curved  channel.  The 
experimenal  study  provided  a  physical  appreciation  and  the 
means  for  theoretical  approximations.  Essentially,  it  was 
found  that  the  jet  behaviour  was  not  significantly  altered 
by  the  presence  of  the  crossflow  except  that  the  jet  was 
deflected  from  a  straight  line  path. 

To  theoretically  predict  the  deflection  of  the  jet  a 
number  of  simpler  problems  were  considered  as  well  as  the 
curved  channel  problem.  The  procedure,  however,  was  the 
same.  A  moment  of  momentum  equation  was  derived  and 
simplified  and  was  found  in  all  cases  to  predict  the 
deflection  given  the  variation  of  the  other  scales  and 
without  the  need  of  additional  empirical  constants. 
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The  theoretical  predictions  were  compared  with  various 


observations  and  were  found  to  be  va 1 i 
portion  of  the  deflected  jet  up  to  a  1 
the  relative  magnitudes  of  the  initial 
longitudinal  velocity,  and  ambient  lat 


for  the  initial 
mi t  which  varies  with 
jet  velocity,  ambient 
;r a  1  veloci  ty . 
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0  .  1  I  nt  roduct i on 

0.1.1  Purpose 

The  effect  of  lateral  currents  on  turbulent  jets  and 
plumes  is  a  complicated  three  dimensional  problem  in  the 
field  of  enviromental  fluid  mechanics.  Much  work  has  been 
done  on  this  problem,  mostly  of  an  experimental  nature,  but 
general  results  are  few.  It  is  the  objective  of  the  present 
study  to  investigate  the  effect  of  weak  lateral  currents  on 
turbulent  jets.  A  specific  practical  problem  to  be 
considered  is  a  circular  turbulent  wall  jet  discharging  into 
a  curved  co-flowing  stream. 

0.1.2  Pr act i ca 1  Si qni f i cance 

Accurate  prediction  of  the  enviromental  effects  of 
industrial  and  domestic  waste  disposal  is  becoming  ever  more 
important  in  the  world  today.  The  first  step  in  predicting 
these  effects  is  normally  to  calculate  the  pollutant 
concentration  field.  The  tools  available  for  this 
calculation  are  jet  theory,  diffusion  theory,  and  dispersion 
theory.  Each  has  it's  own  range  of  applicability  and 
limitations.  In  some  problems  all  three  approaches  may  be 
required  to  cover  the  appropriate  range  of  interest.  The 
present  study  is  concerned  with  the  problems  in  the  domain 
of  jet  theory. 
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Jet  theory  is  most  useful  near  the  pollutant  outlet, 
where  concentrations  are  highest.  The  excess  momentum  of  the 
jet  causes  entrainment  of  the  ambient  fluid  resulting  in 
rapid  mixing  and  dilution.  Bouyancy  effects  may  or  may  not 
be  important.  The  ambient  fluid  body  may  be  a  stream,  a 
river,  lakes  of  various  sizes,  the  ocean,  or  the  atmosphere. 
Rarely,  however,  is  the  ambient  fluid  at  rest.  Thus,  to 
accurately  predict  the  pollutant  concentration  field,  the 
effect  of  the  ambient  motion  on  the  jet  and  pollutant  mixing 
must  be  known.  Essentially,  the  calculation  of  pollutant 
mixing  by  a  jet  is  accomplished  by  calculation  of  the 
velocity  field.  Using  Reynold's  hypothesis  modified  to  fit 
the  experimental  data,  the  concentration  field  may  be 
calculated  from  the  velocity  field  through  the  use  of  a 
pollutant  conservation  equation. 

The  present  study  is  concerned  with  the  discharge  of 
jets  into  a  moving  ambient  fluid  with  a  small  lateral 
component  (relative  to  the  jet).  Some  examples  of  practical 
situations  are:  a  jet  discharging  into  a  stream  at  a  small 
angle  to  the  flow,  a  jet  discharging  into  a  curved  channel 
flow,  a  bouyant  plume  discharged  into  an  ocean  current,  and 
a  heated  surface  discharge  onto  a  lake  with  circulation. 

Although  the  results  of  this  study  may  be  applied  to 
any  of  these  cases,  the  emphasis  will  be  on  the  case  of  a 
discharge  into  a  curved  channel.  As  most  natural  channels 


■ 


3 


are  curved  to  some  extent,  the  resulting  behavior  of  a 
pollutant  discharge  is  of  considerable  practical  importance. 

0.1.3  Sketch  of  Proposed  Work 

The  effect  of  lateral  currents  on  turbulent  jets  was 
investigated  by  the  following  procedure.  First,  an 
experimental  procedure  was  devised  to  obtain  qualitative 
information  useful  in  problem  simplification.  The 
arrangement  used  was  that  of  a  circular  wall  jet  injected 
longitudinally  into  a  curved  channel  flow.  This  setup 
simulated  to  some  extent  the  practical  situation  of  a 
pollutant  diffuser  in  a  natural  channel.  Velocity  profiles, 
both  vertical  and  lateral,  were  taken  at  a  number  of 
downstream  locations  for  a  range  of  initial  velocity  ratios. 
The  position  of  the  jet  in  the  channel  cross  -  sect i on  was 
also  varied. 

The  results  of  the  experiments  were  analysed  with 
respect  to  similarity  of  profiles  and  downstream  scale 
variation.  These  results  were  then  compared  to  the  results 
of  a  similar  study  in  a  straight  channel.  Qualitatively,  the 
effect  of  the  crossflow  was  observed  and  the  important 
parameters  noted. 

Using  the  simplifications  afforded  by  the  experimental 
results,  a  theoretical  analysis  was  made  of  a  similar  but 
more  basic  situation.  The  essential  assumption  justified  by 
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the  experimental  results  was  that  the  cross-sectional 
profiles  of  the  downstream  component  of  excess  velocity 
remained  similar.  With  this  assumption,  integral  momentum, 
energy,  and  moment  of  momentum  equations  were  derived.  These 
equations  provided  the  means  to  calculate  the  the  variation 
of  the  width  scale,  velocity  scale,  and  deflection  of  the 
jet  in  the  downstream  direction.  The  latter  parameter  and 
terms  involving  the  crossflow  were  found  to  occur  only  in 
the  moment  of  momentum  equation  thus  allowing  the  first  two 
equations  to  be  solved  separately  and  the  results 
substituted  into  the  third  to  explicitly  calculate  the 
def 1 ect i on . 

Various  simplified  problems  were  studied  and  the 
results  compared  to  experimental  data  where  available. 
Finally,  the  theory  was  applied  to  the  experimental  work  and 
the  results  compared. 


0 . 2  Review  of  Existing  Work 
0.2.1  Jets  i n  Crossf low 

In  recent  years  the  subject  of  jets  in  crossflows  has 
recieved  considerable  attention  from  various  researchers. 

The  earlier  studies  were  primarily  of  an  experimental  nature 
including  descriptive  observation  and  empirical  predictive 
formulae.  Notable  studies  include  Pratte  and  Baines  (1967) 
Keffer  and  Baines  (1963)  and  Gordier  (1959). 
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Based  on  these  observations  three  separate  regions  of 
the  deflected  jet  may  be  distinguished.  The  first,  the  near 
field,  is  characterized  by  high  jet  velocity  to  crossflow 
velocity  ratio,  and  small  deflections.  The  second  region, 
the  curvilinear  region,  is  char acter i zed  by  large  ambient 
fluid  entrainment,  rapid  jet  velocity  decay,  and  large 
deflection.  The  characteristic  'Kidney'  shaped  velocity 
contours  develop  in  this  region.  The  third  region,  the 
vortex  region,  is  characterized  by  a  pair  of  counterrotating 
vortices  generated  in  the  wake  of  the  deflecting  jet.  By 
this  time  the  jet  axis  is  nearly  parallel  to  the  ambient 
fluid  direction.  In  any  analysis  the  special  features  of 
each  region  must  be  maintained  in  consideration. 

With  reference  to  the  present  study,  attention  is 
focused  only  on  the  initial  near  field  region.  Even  then, 
only  the  special  situation  wherein  this  region  is  much 
longer  than  the  potential  core  length  is  considered. 

The  early  studies  provided  empirical  expressions  for 
the  deflected  jet  axis  usually  in  the  form  of  a  power  law 
ie : 


r  >um 

(  bo  j 


( 1 ) 


where  x  is  longitudinal  distance,  d  is  deflection  ,  s  is  the 
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initial  jet  to  crossflow  velocity  ratio,  b0is  nozzle  radius, 
and  C,  n,  and  m  are  empirical  constants.  Often,  the 
expression  was  inverted  to  give  the  penetration  as  a 
function  of  downstream  distance.  In  some  cases  the  exponent 
and  constant  were  evaluated  separatly  for  the  different 
regions  but  usually  one  formula  was  fitted  to  the  entire 
jet.  A  review  of  these  formulae  appears  in  Rajaratnam 
(1976)  . 

Early  analytical  approaches  to  this  problem,  such  as 
Abramovich  (1963)  and  Crowe  and  Reisbieter  (1967)  were  based 
on  the  principle  that  the  deflection  of  the  jet  is  caused  by 
a  pressure  differential  analogous  to  the  drag  force  on  a 
cylinder.  It  was  found  however,  that  this  mechanism  could 
not  account  for  the  deflection  unless  unrealistically  high 
values  of  the  drag  coefficent  were  chosen. 

A  later  study,  Platten  and  Keffer  (1968),  introduced 
the  mechanism  of  entrainment  as  the  cause  of  jet  deflection. 
Essentially,  the  jet  obtains  momentum  in  the  lateral 
direction  by  entrainment  of  the  ambient  fluid.  This 
entrainment  was  found  to  be  larger  than  for  free  jets  but 
the  form  of  the  entrainment  relation  and  coefficents  proved 
difficult  to  evaluate. 

The  most  recent  analytical  studies  by  Chan,  Lin,  and 
Kennedy  (1976),  Wright  (1977)  and  Makibata  and  IV! i y a i  (1979) 
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all  recognize  the  three  separate  regions  and  use  both  drag 
and  entrainment  concepts  to  calculate  the  deflection  and 
other  parameters.  The  Chan,  Lin,  and  Kennedy  study,  in 
particular ,  uses  integrated  equations  of  motion  and 
similarity  relations  in  the  calculation.  In  regard  to  the 
near  field  region,  however,  the  deflection  is  taken  simply 
as  a  consequence  of  lateral  convection  and  the  equation  for 
the  deflection  is: 


dd 

dx 


V_ 

u 


(2 


m 


Generally,  these  studies  were  primarily  concerned  with 
jets  discharging  perpendicular  to  the  ambient  fluid 
direction  although  some  observations  on  oblique  jets  were 
also  presented.  Bouyancy  effects  were  also  investigated, 
particularly  by  Wright  (1977).  No  study,  however,  included 
wa 1 1  effects . 

0.2.2  T urbu lent  Jets 

The  mechanics  of  turbulent  jets  is  a  well  studied 
problem  and  all  the  essential  information  used  herein  may  be 
easily  obtained  in  Rajaratnam  (1976).  Of  particular  interest 
is  the  problem  of  compound  jets.  The  results  of  an  analysis 
by  Pande  and  Rajaratnam  l  1979)  are  used  in  simplified  form 
in  the  present  study.  Circular  compound  wall  jets  were 
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studied  experimentally  by  Stalker  (1979). 

0.2.3  Curved  Channe 1  F low 

Since  Thomson  (1876)  noted  the  presence  of  a  spiral 
current  in  the  flow  around  bends  of  streams,  much 
investigation  has  been  made  into  the  problem.  The  most 
notable  studies  were  Einstien  and  Harder  (1956),  Rozovski 
(1957),  Yen(1965),  and  Engelund  (1974).  Due  to  the 
complexity  of  the  problem,  these  investigators  were 
concerned  mainly  with  the  special  case  of  developed  or 
ax i symmetr i c  flow.  They  did  however,  obtain  useful 
information  about  the  magnitude  and  shape  of  the  lateral 
velocity  distribution.  In  particular,  the  approximate 
relation  due  to  Engelund  (1974): 

V_  „  7  h  (3) 

U  R 

where  h  is  the  depth  of  flow,  R  is  the  radius  of  curvature, 

U  is  mean  downstream  velocity,  and  V  is  lateral  velocity 
gives  an  estimate  of  the  magnitude  of  the  crossflow  near  the 
bed.  Over  the  depth  of  the  jet  the  average  velocity  ratio  is 
given  by  changing  the  constant  to  about  three. 


1.  EXPERIMENTAL  INVESTIGATION 


1 . 1  Exper iments 

1.1.1  Apparatus  and  Equipment 

1  .  1  .  1  .  1  F 1 ume 

The  flume  used  was  curved  in  an  'S'  shape  with  two 
consecutive  180  degree  bends  of  9.0  feet  centerline  radius. 
The  test  section  was  located  at  90  degrees  into  the  second 
bend.  A  plan  of  the  flume  is  shown  in  figure  (1)  and  a  view 
of  the  test  location  is  shown  in  plate  (1).  The  flume  was 
located  at  the  Ellerslie  River  Engineering  research  station 
of  the  University  of  Alberta. 

The  flume  is  equiped  with  it's  own  sump  and  pump  and 
has  an  adjustable  tailgate  to  regulate  flow  depth.  The  table 
slope  of  the  flume  is  adjustable  but  was  mantained  level  for 
the  experiments.  The  bed  was  smooth  aluminum  and  was  marked 
with  a  curvilinear  grid  for  qualitative  reference. 

1  .  1  . 1  . 2  Traverse 

The  traverse  used  provided  electrically  powered  motion 
in  the  lateral  and  vertical  directions  but  required  manual 
adjustment  in  the  longitudinal  direction.  The  lateral 
position  was  measured  electronicly  by  means  of  a 
potentiometer  to  the  nearest  one  hundredth  of  a  foot  while 
the  vertical  position  was  measured  visually  using  a  staff 
gauge  with  a  vernier  to  the  nearest  one  thousandth  of  a 
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Figure  1  -  Plan  of  Curved  Flume 


Plate  1  -  Test  Location 
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foot.  The  downstream  locations  were  set  at  0.2  foot 
intervals  along  lines  tangent  to  the  curve  at  the  jet 
location.  Plate  (2)  shows  the  traverse. 

1  . 1 . 1 . 3  Jet  Arrangement 

The  jet  arrangement  involved  a  constant  head  device, 
large  diameter  tubing,  and  a  nozzle.  The  constant  head 
device  was  mounted  on  an  adjustable  stand  and  fed  by  a  hose 
directly  from  the  city  water  supply.  The  tubing  was  3/4  inch 
inside  diameter  plastic  ' Tygon  Tubing'  and  was  about  five 
feet  in  length.  The  nozzle  was  a  steamlined  tube  bent  at  a 
90  degree  angle  at  the  bottom  to  the  exit.  It  was  introduced 
into  the  flow  through  the  water  surface  to  the  bed.  This 
resulted  in  some  disturbance  to  the  flow  but  was  necessary 
because  of  the  requirement  that  the  location  of  the  jet  be 
changed  easily.  The  inside  diameter  of  the  tube  portion  was 
3/8  inch  and  the  outlet  diameter  was  1/4  inch. 

The  elevation  of  the  water  in  the  head  tank  was  not 
measured  because  it  was  found  that  the  losses  in  the  tubing 
and  nozzle  were  very  significant  and  therefore  the  only 
reliable  measurement  of  initial  velocity  was  an  actual 
velocity  reading  in  the  potential  core  of  the  jet. 

Two  flaws  of  this  setup  became  apparent  during  and 
after  the  course  of  experiments.  The  first  was  that  the  jet 
flow  was  at  city  water  temperature  while  the  channel  flow 
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Plate  2  Traverse  for  Velocity  Measurements 
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tended  to  be  at  room  temperature,  resulting  in  possible 
bouyancy  effects.  The  second  problem  was  that  there  was  no 
provision  to  accurately  align  the  jet  with  the  coordinate 
system.  As  a  consequence  the  measured  deflections  proved 
difficult  to  correlate. 

1 ■ 1 ■ 1 . 4  Measurement  System 

Velocity  measurements  were  made  by  means  of  a  standard 
Prandtl  tube  of  1/8  inch  outside  diameter.  The  pressure 
difference  was  measured  with  a  sloping  manometer  set  at  0, 
60,  and  75  degrees  from  the  vertical  depending  on  the 
magnitude  of  the  velocity.  At  the  flatest  setting,  pressure 
differentials  could  be  measured  to  the  nearest  one 
thousandth  of  a  foot  of  water. 

Water  surface  slope  was  measured  by  means  of  two  static 
head  screwdriver  type  probes  placed  at  45  degrees  upstream 
and  downstream  of  the  jet  location.  The  depth  at  these 
locations  was  indicated  on  the  same  sloping  manometer  as  the 
velocity  readings. 

A  yaw  probe  was  also  considered  for  use  but  the 
additional  information  about  the  lateral  velocities  was 
found  to  be  unreliable  due  to  the  width  of  the  probe 
compared  to  the  lateral  velocity  gradient. 
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1.1.2  Procedure 

1 . 1 . 2  .  1  Pre 1 i mi  nary  Runs 

A  series  of  preliminary  runs  were  performed  in  order  to 
qualitatively  examine  the  problem.  These  runs  used  the 
equipment  described  above  but  no  velocity  measurements  were 
made  other  than  to  obtain  an  initial  velocity  ratio. 

Instead,  the  injected  water  was  dyed  and  the  dispersion  of 
the  jet  was  visually  studied.  A  typical  plot  of  a  visually 
time  averaged  jet  profile  is  shown  in  figure  (2).  A 
photograph  of  the  dyed  jet  is  shown  in  plate  (3).  Note  that 
the  plot  is  referenced  to  the  arc  of  the  curve  determined  by 
the  location  of  the  jet.  Various  parameters  such  as  velocity 
ratio,  channel  Reynolds  and  Froude  numbers,  and  lateral 
position  of  the  jet  were  varied.  The  parameters  selected  for 
further  study  were  the  velocity  ratio  and  jet  position. 

1 . 1 . 2 . 2  Actua 1  Runs 

After  completion  of  the  preliminary  runs  a  series  of 
six  more  detailed  runs  were  performed.  Table  (1)  indicates 
the  appropriate  parameters  for  each  experiment.  Figure  (3) 
is  the  definition  sketch  for  this  data.  The  object  of  these 
runs  was  to  obtain  qualitative  and  quantitative  information 
about  the  velocity  distribution  of  the  jet. 

Each  experiment  involved  taking  a  lateral  and  vertical 
velocity  profile  at  three  or  four  downstream  stations.  The 
only  component  of  velocity  measured  was  the  '  u'  component; 
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Figure  2  -  Plot  of  the  Boundaries  of  Dyed  Jet  (Plan  View) 


17 


Plate  3  -  Plan  View  of  Dyed  Jet 
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Figure  3.  Experiment  Definition  Sketch. 
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that  is,  the  component  parallel  to  the  initial  jet  axis.  The 
lateral  profile  was  taken  first  with  the  Prandtl  tube 
slightly  above  the  flume  bed.  The  results  of  the  lateral 
profile  were  plotted  as  they  were  taken  and  the  position  of 
maximum  velocity  noted.  About  ten  points  were  found  to  be 
sufficent  to  describe  the  lateral  profile. 

The  vertical  profile  was  taken  at  the  location  of 
maximum  velocity  as  indicated  by  the  lateral  profile.  The 
first  reading  was  taken  with  the  tube  resting  on  the  bed  and 
subsequent  readings  were  taken  at  spacings  upward  to  allow 
the  profile  to  be  adequately  defined  with  approx i amet 1 y  ten 
points.  The  elevation  of  the  lateral  profile  was  checked  to 
insure  that  it  was  reasonably  close  to  the  elevation  of 
maximum  velocity. 

To  insure  that  air  bubbles  did  not  interfere  with  the 
readings,  the  system  was  flushed  before  each  profile. 
Readings  were  then  taken  in  order  of  decreasing  velocity  to 
prevent  the  air  rich  ambient  water  from  entering  the  system. 
For  the  lateral  profiles  this  involved  alternating  the  pitot 
tube  from  side  to  side  of  the  centerline  and  further  away. 

Figures  ( 4 )  - ( 15)  show  the  vertical  and  lateral  velocity 
profiles  for  all  sections  of  runs  one  to  six. 
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Figure  4  Run  1  Lateral  Profiles  of  Longitudinal  Velocity 
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Figure  5  -  Run  1  Vertical  Profiles  of  Longitudinal  Velocity 
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Figure  6  -  Run  2  Lateral  Profiles  of  Longitudinal  Velocity 
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Figure  7  -  Run  2  Vertical  Profiles  of  Longitudinal  Velocity 
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Figure  8  -  Run  3  Lateral  Profiles  of  Longitudinal  Velocity 
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Figure  9  -  Run  3  Vertical  Profiles  of  Longitudinal  Velocity 
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Figure  10  -  Run  4  Lateral  Profiles  of  Longitudinal  Velocity 
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Figure  11  -  Run  4  Vertical  Profiles  of  Longitudinal  Velocity 
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Figure  12  -  Run  5  Lateral  Profiles  of  Longitudinal  Velocity 


' 


. 

CD  CD  CD  CD 
CM  CD  OO 

mama 

OCDOO 
II  II  II  II 
X  X  X  X 


B0<< 


CJ 

LU 

CO 

\ 

♦ 

(— 

u_ 


3 


Figure  13  -  Run  5  Vertical  Profiles  of  Longitudinal  Velocity 
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Figure  14  -  Run  6  Lateral  Profiles  of  Longitudinal  Velocity 
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Figure  15  -  Run  6  Vertical  Profiles  of  Longitudinal  Velocity 
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1 . 2  Ana  lysis  of  Exper iments 

1.2.1  Methodology 

To  reduce  the  experimental  data  to  a  managable  and 
useful  form,  similarity  analysis  was  applied.  Practicaly, 
this  involved  plotting  each  profile,  fitting  a  smooth  curve 
through  the  data  points  and  then  extracting  the  pertinent 
scales.  The  scales  needed  were  maximum  velocity,  distance  of 
maximum  velocity  from  origin,  and  distance  from  point  of 
maximimum  velocity  to  the  point  where  the  velocity  was  one 
half  maximum.  In  the  case  of  lateral  profiles  the  latter 
scale  was  measured  in  both  directions  and  an  average  taken. 

Using  these  scales,  which  were  unique  to  any  one 
profile,  the  similarity  profiles  of  each  run  were  plotted 
together  and  were  indeed  found  to  be  similar.  Figures  (16)  - 
(27)  show  these  similarity  plots  for  all  the  runs.  The  solid 
line  is  the  often  used  exponential  curve: 
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Figure  16  -  Lateral  Similarity  Profile  of  Longitudinal 
Ve loci ty  Run  1 
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Figure  17  -  Vertical  Similarity  Profile  of  Longitudinal 
Ve loci ty  Run  1 
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Figure  18  -  Lateral  Similarity  Profile  of  Longitudinal 
Ve loci ty  Run  2 
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Figure  19  -  Vertical  Similarity  Profile  of  Lonqitudinal 
Ve loci ty  Run  2 
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Figure  21  -  Vertical  Similarity  Profile  of  Longitudinal 
Ve loci ty  Run  3 
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Figure  22  -  Lateral  Similarity  Profile  of  Longitudinal 
Veloci ty  Run  4 
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Figure  23  -  Vertical  Similarity  Profile  of  Longitudinal 
Ve loci ty  Run  4 
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Figure  24  -  Lateral  Similarity  Profile  of  Longitudinal 
Veloci ty  Run  5 
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Figure  25  -  Vertical  Similarity  Profile  of  Longitudinal 
Ve loci ty  Run  5 
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Figure  26  -  Lateral  Similarity  Profile  of  Longitudinal 
Ve loci ty  Run  6 
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Figure  27  -  Vertical  Similarity  Profile  of  Longitudinal 
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1.2.2  Sea  1 e  Ana  lysis 

The  velocity  data,  having  been  shown  to  have  similar 
profiles,  were  next  analysed  with  respect  to  the  variation 
of  their  scales.  Figures  (28)  -  (32)  show  graphically  the 
downstream  variation  of  each  of  the  scales  for  all  of  the 
runs.  In  order  to  correlate  these  scale  variations  reference 
was  made  to  the  compound  jet  theory  briefly  discussed  above. 
Using  U  as  a  velocity  scale  and  the  momentum  thickness  -  to 
nondimens i ona 1 i se  the  profile  scales,  figures  (33)  -  (37) 
were  plotted  and  show  reasonably  good  corelation.  The 
results  of  Sta Iker ( 1 979 )  are  shown  for  comparison.  The 
velocity  correlation  agrees  very  well  while  the  length  scale 
predictions  do  not.  This  is  believed  to  have  been  caused  by 
the  small  density  difference  between  the  jet  and  ambient 
flow  which  resulted  in  a  flatter,  wider  jet  profile  in  the 
present  study. 
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Figure  28  -  Variation  of  Ma 
Jets  in  Curved  Channel  Flow 


RUN 


48 


— '  CM  00  ^  LO  CO 


=3  ZD  ZD  ZD  ZD 
O'  Q£  C3C  QZ.  O' 


©  <  + 


B  +  <  04 


4  E  €)  X 


©# 


B 


£_ 

(0 

o 

■ — 

LO 

D 

O 

m 

c_ 

o 

■r— 

o 

1 

1 

<D 

LO 

■ — 

CM 

CO 

CO 

o 

m 

OO 

o 

-C 

4-» 

D> 

c 

0) 

o 

_J 

CD 

LO 

r— 

m 

CO  2 

o 

0  0 

■  r-  r— 

*->  U_ 

£_ 

in 
oo  *x* 


O 

LO 

CSJ 


T - T 

2£0*  0 


LO 

CM 


CD 

CD 

CD 


frZO*  0  910*0 

rid)  zq 


a)  — 

>  (D 

c 

M-  C 
O  to 
_c 
c  o 
o 

—  T3 
+->  0) 
CO  > 

•r-  [_ 

£-  D 
CO  O 


CD  U) 
CM  -*-< 
<D 

<1)  ~0 
C_ 

D  — 
O)— 

•r-  (0 


OtO'O 


- 1 — 

800*0 


000*0 


49 


•f- 

OSI-O 


CM  CO  -^LO  CD 

Z22ZZ2 
ZD  ID  ZD  ZD  ZD  ZD 

occnoccuctiaL 


BQ< +  XO 

- * - 1 - 1 - 1 - -I - -i 

©  ^  +  EM 


□  t>  X  +  ◄ 


E*>«  e  x 


^ - 1 - 1 - ( - 1 - 1 - r 

021  *  0  060*0  090*0 

VLB 


LO 

1 - -I - -f-oo 

m 

o 


*€©-E 


CD 

£_ 

to 

LO 

' — 

-p- 

15 

9 

O 

CD 

c_ 

LO 

o 

1 

0) 

CM 

■ — 

-  CO 

(0 

9 

o 

CD 

to 

CD 

JZ 

4-J 

05 

c 

<D 

CD 

_j 

-  LO 

9 

—  3 

CD 

to  0 

♦ 

£_  «— 
<D  U_ 

1 - 

Ll 

CO  — 

—i  0) 

lO 

c 

X 

M-  C 

-CD 

0  co 

9 

_c 

CD 

C  CD 

CD 

0 

"O 

-t->  (1) 

to  > 

•r-  C_ 

£-  D 

LO 

CO  CD 

-  CM 

> 

9 

C 

O 

1  'I— 

LO 

re  30 
Jets 

<SJ 

=5  — 

“  vH 

o>— 

9 

(0 

O 

Ll.  32 

- 1 - 1 - 

060*0  000*0 


CD 

CD 

CD 

9 

O 


50 


CM  CO  ''t  LO  CO 


^rDZ3IOZ3ZD 
od  cm  cm  a:  oc  on 


BQ<+XO 


010*0 


000*0 


31  -  Variation  of  Boundary  Layer  Thickness  -  Circular 
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Figure  32  -  Jet  Deflection 
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Jets  in  Curved  Channel  Flow 
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2.  THEORETICAL  INVESTIGATION 

2 . 1  Theoret i ca 1  Analysis 

2.1.1  Pre 1 i mi  nary  Remarks 

The  experimental  work  above  showed  that  the  jet 
parameters  and  similarity  profiles  were  not  significantly 
altered  by  the  presence  of  the  crossflow.  The  jet,  however, 
was  deflected  to  follow  a  curved  path.  In  this  chapter  a 
theoretical  analysis  based  on  an  assumption  of  a  weak 
crossflow  will  be  presented. 

To  start  with,  for  simplicity,  the  case  of  a  plane, 
turbulent  jet  in  a  uniform  crossflow  will  be  studied.  The 
technique  developed  will  then  be  applied  to  the  more 
complicated  problems  of  a  circular  compound  jet  as  well  as 
to  the  circular  wall  jet  in  a  curved  stream. 

2.1.2  P 1 ane  Jet  i n  Un i form  Crossf 1 ow 

Refering  to  figure  (38),  the  definition  sketch  for  a 
plane  jet  in  a  crossflow,  the  situation  is  that  of  a  plane 
jet  of  initial  velocity  U0and  width  2b0issuing  into  an 
ambient  fluid  with  a  velocity  V  in  a  direction  perpendicular 
to  the  jet  axis.  The  parameters  which  characterize  the 
velocity  field  are:  the  maximum  '  x'  component  of  velocity  at 

any  section,  u  ,  the  half  velocity  width  of  the  jet  b  ,  and 

m 

the  deflection  of  the  jet,  d,  defined  to  be  the  distance 
from  the  outlet  axis  to  the  location  of  velocity  maximum. 


. 
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Figure  38.  Definition  Sketch  of  Plane  Jet  in  Crossflow. 
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Assuming  V  is  of  the  order  of  magnitude  of  v'  and  the 
pressure  gradient  terms  are  negligable,  the  equations  of 
motion  for  this  problem  reduce  to  (see  the  appendix  for 
reduct  ion ) : 


u~  +  v  •  ly.  +  v—  =  A  <LL 

9x  9y  9y  p  9y 


(5) 


iH.  +  9v'=  n  (6) 

9x  9y 


Where : 


v'=  v  -  V 


(7  ) 


The  prime  on  v'  will  be  dropped  from  this  point  on  for 
convi eni ence .  Integral  relationships  will  now  be  developed 
to  predict  um ,  b  ,  and  d. 


2  .  1  .  2  .  1  Integra  1  Momentum  Equat ion 

Integrate  equation  (5)  with  respect  to  y  from  minus 
infinity  to  plus  infinity  to  obtain: 


oo 


9  u  v 
9y 


dy  + 


9u  V 

9y 


dy  =  -p^ 


(81 


-QC 


— OO 


— oo 


—00 


' 
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since  u  and  t  equal  zero  at  plus  and  minus  infinity,  this 
equation  reduces  to: 


d_ 

dx 


u2  dy  =  0 


(9) 


This  equation  states  that  total  axial  momentum  flux  is 
conserved . 

2  .  1  .  2 . 2  Integral  Energy  Eguat i on 

Multiply  equation  (5)  through  by  u  and  use  (6)  to 
rearrange  it  into: 


9u3  9u2  v  +  9u2V  _  2  9t 
9x  9y  9y  p  9y 


10) 


Again,  integrate  with  respect  to  y  as  before 


f  3U^y  +  [  + 

J  9x  • y  J  9y 


9u2V^„  _  2  f 


9t 


dy  =  “  j  u  sy 


9y  p 


1 1 


Simplifying,  the  integral  energy  equation  is  obtained. 


oo 


dx 


u3dy 


P  J 


u!ydy 


(  12) 


This  equation  states  that  the  rate  of  change  of  total 


' 
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Kinetic  energy  flux  in  the  jet  is  equal  to  the  dissipation 
caused  by  the  shear  stress. 


2 . 1 . 2 . 3  Comment 

The  transverse  velocity  V  has  disappeared  from  both  the 
momentum  and  energy  equations  which  would  seem  to  imply  that 
it  has  no  effect  on  the  axial  momentum  and  energy  balances. 


2 . 1 . 2 . 4  Integral  Moment  of  Momentum  Equat ion 

Multiply  equation  (5)  through  by  y  and  integrate  as 
above : 


3y  u- 


3x 


-dy  + 


3uv  , 

yVy  + 


3uV ,  1 

yVy=  P 


3t  , 

y5ydy 


(13 


-00 


Integrating  the  second,  third,  and  fourth  terms  by  parts 


rl  f  f 

.  y  u2  dy  -  j  u  v  dy  - 


u  V  dy  =  -M  t  dy 


(  14  ) 


Since  u  goes  to  zero  exponentially  the  terms  evaluated  at 
infinity  are  all  equal  to  zero.  The  remaining  term  involving 
shear  stress  is  also  zero  by  antisymmetry.  This  leaves: 


oo 


d_ 

dx 

—00 


y  u2dy  - 


y  v  dy  - 


u  V  dy  =  0 


(  15) 


-00 


-00 


. 
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Consider  the  second  and  third  terms.  It  may  be  noted  that  in 
the  second  term  v  will  be  small  where  u  is  largest  and  also 
that  v  wi 1 1  tend  to  have  opposite  signs  on  opposite  sides  of 
the  jet.  Thus  the  integral  will  tend  to  be  small.  If  V  is  of 
the  same  order  of  magnitude  as  v  then  the  third  term  will 
clearly  be  much  larger  and  the  second  may  be  neglected. 
(Note:  In  an  undeflected  jet  v  is  antisymmetric  and  u  is 
symmetric  resulting  in  the  integral  going  to  zero 
identically.  In  the  weakly  deflected  case  considered  here  u 
remains  symmetric  but  v  is  no  longer  antisymmetric.  The 
deviation,  however,  is  small.)  Thus  the  following  is 
obtained . 


co 


CO 


(  16  ) 


This  equation  may  be  interpreted  as  stating  that  the  rate  of 
change  of  the  total  flux  of  moment  of  momentum  equals  the 
convection  of  u  velocity  in  the  y  direction. 

2 . 1 . 2 . 5  Simi lar i ty  Hypothesis 

Following  the  experimental  analysis  and  preliminary 
remarks  the  profile  of  u  velocity  in  the  y  direction  is 
assumed  to  be  similar,  or: 


(17) 


. 

■ 
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Also  assume: 


T 


Mn) 


(  18) 


Note  that: 


y  =  nb  +  d 


(  19) 


dy  =  b  dri 


(20) 


Substitute  equations  (17),  (18),  (19),  and  (20)  into 
equations  (9),  (12),  and  (16)  to  obtain: 


d  u  2 
dx  m 


00 

b  [  f ! 2  din  =  0 
J 

—00 


(21  ) 


,  „3b  f:3dn 
dx  m  J 


2  Um3 

m 


fi  f 2 1 dn 


(22) 


-oo 
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oo 


oo 


d_ 

dx 


(n  b  +  d)  u  2 f 2 2 b  dn  = 
m 

—oo 


f  b  dri 


(23) 


Simplify  to  get : 


d_ 

dx 


(u  2b)  =  0 
m 


(24) 


d_ 

dx 


oo 


u  3b 
m 


fi3dn  =  2  u  3 
m 


fi 


-oo 


f  2  '  dri 


(25) 


d 

dx 


(u  2b  d) 
m  ‘ 


Li 

Iz 


V  b  u 

m 


(26) 


It  may  be  noted  that  equations  (24)  and  (25)  are  identical 
to  the  corresponding  equations  for  an  undeflected  jet.  Since 
these  equations  do  not  contain  d  or  V  they  may  be  solved  to 
give  umand  b  .  The  solution  will  be  identical  to  the  case  of 
an  undeflected  jet  except  for  a  possible  difference  in  the 
empirical  constant.  Assuming  for  now  that  the  constants  will 
indeed  be  the  same,  Ra jaratnam( 1 976 )  gives  the  solutions  as: 


3.50 


(27) 


/  x  /  b0 


' 


or  tsup* 
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b  =  0.097  x 


(28) 


These  expressions  may  now  be  substituted  into  (26)  to 
calculate  d.  Letting: 


Ii  = 


fi  dn 


29) 


I2 


—oo 


(30) 


and  substituting  (27)  and  (28)  into  (26) 


TL  l3-:-5-  u.o.)2  b_o..X°.-0_97 .  x ).  d  =  h-  u0  V  (0.097  x)  [  — --5— - 

QX  X  1 9  /  /  l 

(/x  /  b0 J 


(31  ) 


Rearrangi ng : 


TP-  =  0.286 

dx 


i l  )L  n, 

I  2  U  o  /  b0 


(32) 


Solving  the  differential  equation  with  the  initial 
cond i t i on : 


9 


* 


66 


d  =  0  at  x  =  0 


(33) 


we  obtain: 


0.19  TT-  4-1- 


V  Ii 


r  \ 

x 


Uo  I 


2  !^0j 


3  /  2 


(34) 


Evaluating  constants  using  the  convienient  exponential  form 
for  f : 


f  (ti) 


-0.693  n2 
e 


(35) 


g i ves : 


Ij  =  1.063 


(36) 


I2  =  0.751 


(37) 


Finally,  substituting  for  these  constants  we  obtain: 


^3/2 


(38) 
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Equation  (38)  gives  the  deflection  of  a  simple  plane  jet  as 
a  function  of  downsteam  distance  and  initial  jet  parameters. 

2.1.3  Ci rcu 1 ar  Compound  Jet  In  Uni  form  Crossflow 

The  technique  for  analysing  a  circular  compound  jet  in 
a  uniform  crossfow  is  identical  to  that  used  to  analyse  the 
plane  jet.  The  symmetry  conditions,  however,  are  no  longer 
satisfied  and  it  will  be  necessary  to  integrate  in  both 
cross-sectional  directions.  The  more  complicated  boundary 
conditions  at  infinity  will  also  result  in  added  complexity 
in  both  the  analysis  and  final  expression  for  the 
deflection.  Refering  to  the  definition  sketch,  figure  (39), 
the  additional  parameters  are  the  background  axial  velocity 
U  and  b  is  now  the  nozzle  diameter.  The  equations  of  motion 
for  this  problem  are: 


9u  +  3v  +  8w=0  (39) 

3x  3y  3z 


3y 


3u 

3z 


vi“  =  I 

3y  p 


3t 


zx 

3z 


I 

. 

I 

J 


+ 


(40) 
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x 


Figure  39.  Definition  Sketch  of  Circular  Compound  Jet  in  Crossflow. 
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2 . 1 . 3 . 1  Integral  Momentum  Equa t i on 

Integrate  (40)  with  respect  to  '  z'  from  minus  to  plus 
infinity. 


00 


3 

¥x 


u2dz  + 


f  auvj  + 


J  1? 


8x 


U  u  dz  + 


zx 


9y 


V  dz  =  - 


f  3t 


yx 


9y 


dz 


(42) 


-00 


Now  integrate  (42)  with  respect  to  '  y1  from  minus  infinity 
to  plus  infinity  and  simplify  as  above. 


djf 

dxjj 

— oo 


u  (U  +  u)  dy  dz  =  0 


(43) 


This  equation,  the  integral  momentum  equation,  states  that 
the  exess  momentum  efflux  of  the  jet  is  conserved. 


2 . 1 . 3 . 2  Integral  Energy  Eguat ion 

Multiply  equation  (40)  by  ' u'  and  use  (39)  to  obtain: 


£u3  3u2v  8u2w  +  9U_u2j.  8V  u2^  2u 
8x  +  3y  3z  8x  3y  p 


8t 


xy 


9t. 


9y 


+ 


zx 


9z 


(44) 


Integrating  with  respect  to  '  z'  and  '  y'  as  above  and 
simplifying  the  integral  energy  equation  is  obtained. 


d_ 

dx 


00 

r  r  f 

00 

9t  9t 

u2(U  +  u)  dy  dz  = 

J  J  ^ 

2u 

yx  +  zx 

J  P 

ay  j 

dy  dz 


(45) 


•  r  : 

. 
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This  equation  states  that  the  rate  of  change  of  total  energy 
efflux  equals  the  dissipation  by  shear  stresses. 

Note  once  again  that  the  crossflow  V  does  not  appear  in 

these  equations.  If  we  assume  that  the  width  scales  are  the 

same  in  both  cross-sectional  directions,  then  these  two 

equations  can  be  solved  for  u  and  b.  Note  also  that  these 

m 

equations  are  the  same  as  the  equations  derived  by  Pande  and 
Ra jaratnam( 1 979 )  (though  in  a  different  form). 


2 . 1 . 3 . 3  Integral  Moment  of  Momentum  Equat ion 

Multiply  equation  (40)  by  ' y'  and  integrate  wi th 
respect  to  ' z'  from  minus  infinity  to  plus  infinity. 


00 


8Uuy 

3x 


.  A  f  8Vu  , 

dz  +  j  y  dz 


8t 


=  y 


yx 


ay 


dz 


(46) 


Integrate  again  with  respect  to  ' y'  from  minus  infinity  to 
plus  infinity. 


^  yu  (u  +  U)  dy  dz  + 


8uv  ,  .  ,  f 

y  tc  dy  dz  + 


—  00 


00 


1 

p 


ay 


3Tvx 

y  — ~  dy  dz 


JJ 


avu  ,  , 

y  t —  dy  dz 
J  ay  J 
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Integrating  the  second  term  by  parts: 


0°  , 

1^-  dy  =  [y  u  v]  -  j  u  v  dy 

"  _C°_oo 
0 


(48) 


.* 


'  V 
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Similarily  the  third  term: 


oo 


[y  v  u] 


00 

* 

u  V  dy 

—00 


(49) 


And  the  fourth  term. 


oo 


y 


xy 


9y 


dy 


— oo 


y  t 


yz 


00 


—00 


(50) 


Collecting  the  remaininq  terms: 


00 


d_  f 
dx 

—00 


y  u  (U  +  u)  dy  dz  -  j" 


oo 

f 


y  v  dy  dz  - 


u  V  dy  dz 


J  «/ 

—00 


(51 


€0 


dz 


The  last  term  may  be  set  equal  to  zero  by  antisymmetry.  As 
in  the  plane  jet,  the  second  term,  while  not  identically 
zero,  is  small  compared  to  the  third  term  and  may  be 
neg 1 ected .  Finally: 


d_ 

dx 


J  J 


y  u  (U  +  u)  dy  dz  = 


_  f 


J 


u  V  dy  dz 


-oo 


(52) 


2 . 1 . 3 . 4  Simi 1 ar i ty  Hypothes i s 


Assume : 


. 
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9i 


=  fi(c)  9i(n) 


(53) 


where  fl  and  gl  are  the  similarity  profiles  in  the  z  and  y 
directions.  Note  that  the  width  scale  b  is  assumed  to  be  the 
same  in  both  directions.  Also  assume: 


^  2  -  ^ 2  (  C  )  g  2  (ci  ) 
p  m 


(54) 


T 


P 


ZX 

U 


=  f3U)  93  (p) 


(55) 


Now  substitute  these  similarity  expressions  into  equations 
(43) , (45)  and  ( 52 ) . 


oo 


-CO 


r 


u 


u  b2  f ! gi  !  l  +  -n^  figi 


m 


U 


dri  dc,  =  0 


(56) 


d_ 

dx 


* 

r  r 

u 

2b2 

fi2gi2 

m 

; 

J 

V 

-OO 

u 

+  JH  f 

u  1 


1 

dn  dc 


OO 


(57) 


-oo 


fl  gi  (f 2  92*+  f 3 ' gs )  dp  dc 
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oo 


U  i  2 

-r-  u  b 

dx  m 


fi  gi  d  (U  +  um  fi  gi)  dri  dc  =  V  u  b2 

m 

.00 


f i  gi  do  dc  ( 58 ) 


Equations  (56)  and  (57)  are  found  to  be  identical  to  the 
corresponding  equations  for  an  undeflected  jet.  The 
parameters  umand  b  are  therefore  given  by  the  analysis  of 
Pande  and  Ra jar  a tnam( 1 979 )  as: 

for  I1  <  150  ( 

V 


U__  1  x 

u  ~  k  e 


b 


in 

m 


m 


U 


F?  + 


U 


Fi 


(60) 


where 


(61) 


is  the  momentum  thickness  and  K  is  a  constant  taken  to  be 
equal  to  12.17.  Rewriting  equation  (58) 


d_ 

dx 


(u 


m 


b2d  (U,  U  +  um  I s )1  "  V 


u  b2  I 
m 


(62) 


. 
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where 


00 

•  f 

fi  gi  dn  dc 

J  . 

—00 


(63) 


I5 


(64) 


Substituting  (59)  ,  (60),  (63)  and  (64)  into  (62): 


d  _ x  d 

dx  K  F2  +  ( x/9 )  F7 


K  I5] 

(xTeTj 


v  u  (x/e) 

U  (K  F,  +  (x/e) Fn ) 


(651 


Expand  derivative: 


dd  d  K  (F?  I4  -  Fi  I c) _  V  I4  (x/9)  (66) 

dx  e  (f:  (x/e)  +  k  f2)(i4  (x/e )  +  k  i5)  u  k  i5  +  i4  (x/e) 


Using  equation  (4)  for  both  fi  and  gi  we  may  evaluate  the 
constants  as: 

e-0.693(n2  +  C2) 


—CD 


dn  dc  =  4.533 


(67) 


■ 
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Is 


’  e~0. 693(2) (n2  +  c2) 


dn  dc  =  2.266 


(68) 


From  Pande  and  Rajaratnam  (1979) 


Fi  =  0.7215 


(69  ' 


F2  =  0.3608 


(70) 


therefore : 


(71  ) 


Thus  we  obtain  the  simplified  form: 


dd  _  V  (x/6)  (72) 

dx  "  U  ( (x/0)  +  (K/2) ) 


This  may  be  easily  integrated  to  give: 


d 

e 


v 

u 


x  K  p  f(x/6)  +  (K/2) 
6  "  2  tnl  (K72) 


(73) 


76 

Again,  we  have  a  simple  expression  for  the  calculation  of 
jet  deflection. 


2 • 1 • 3 . 5  Simple  Jet  As  a  Spec i a  1  Case 
Take  equation  (72)  again: 


dd  _  V  (x/6) 

dx  '  U  C(x/G)  +  (K/2) ) 


(74) 


oi  nee 


0  = 


75 


We  have: 


dd  =  _ V  x _ 1 _ 

/■' I  v  - —  -  ■■  ■  ..  —  f  *> 

x_  1  +  K 

r°  /urm  rr  2 

/  U  [u  i 

Setting  U  equal  to  zero: 


r0  U 


U 


Uo.  i 
U  "  ) 


(76) 


dd  _  2  V  x 

dx  K  r0  U0 


(77) 


Solving: 


d_  =  1  V_ 
r o  K  Uo 


fc-]2 

(r0j 


(78) 
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2.1.4  Jet  Pi scharqi nq  On  Bed  Of  Curved  Channel 

The  theory  developed  above  will  now  be  applied  to 
predict  the  jet  deflection  in  a  complicated  practical 
situation.  The  problem  will  be  idealized  as  a  three 
dimensional  compound  wall  jet  with  a  linearly  varying 
crossflow.  In  this  case  no  analytical  solution  exists  that 
the  writer  is  aware  of  for  the  undeflected  case.  However, 
the  experimental  results  of  both  the  present  study  and  that 
of  Sta lker ( 1 979 )  may  be  used  instead  to  provide  the 
necessary  information  about  the  variation  of  width  and 
velocity  scales.  Therefore,  only  the  integral  moment  of 
momentum  equation  will  be  derived  and  used. 


2 . 1 . 4  .  1  Integral  Moment  of  Momentum  Equat i on 

Multiplying  equation  (40)  by  '  y'  and  integrating  with 
respect  to  '  z'  from  z  =  0  to  z  =  00  we  obtain: 


00  00  oo 


Integrate  now  with  respect  to  y  from  y  =  -00  to  y  =  00 


CC  00  00  30 

^  j  j  y  u2dy  dz  +  j|  y  dy  dz  + 

_cc  0  -oo  0 

oc  oo 


00  00 
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00  00 


3^  dy  dz  + 

3X  Jj 
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SuV 


—oo  0 


-  dy  dz 
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p  j  r  zx 

-oo  0 


dy 


Evaluating  the  terms  as  above: 
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OO  oo 


00  00 


y  u  (U  +  u)  dy  dz  - 

J  • 

—oo  0 


u  V  dy  dz  =  0 


(81  ) 


Thus  the  integral  moment  of  momentum  equation  is  identical 
to  that  of  the  previous  section.  In  this  case  however  the 
two  length  scales  may  not  be  taken  as  equal. 


2 . 1 . 4 . 2  Simi 1 ar i ty  Hypothes i s 

Substituting  the  similarity  hypothesis  into  equation 

(81): 


00  oo 

r  f 


-T-  u  b  b 
dx  my  z  . 

u  j  J 

— oo  0 


fi  9i  d  (U  +  fi  gi)  dn  dc 


(82) 


or  : 


=  V  u  b  b  |[  fi  gi  dn  dc 
m  z  y  jj 

—oo  0 


4—  fu  b  b  d  (U  I6+  u  I 7 ) 1  -  V  u  b  b  I6 

dx^myz  m  j  mzy 
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where 


Is 


CO  00 


JJ 


_co  0 


fl 


gi  dn  dc 


(84) 


I  7 


CO  oo 

f!2f22dn  dc 


_ oo  0 


(85) 
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To  solve  for  d  we  will  assume: 


(86) 


(87) 


U 


(88) 


where  the  K  '  s  are  to  be  determined  from  experiment. 
Nondimens iona 1 ize  variables  by  letting 


(89) 


Substitute  the  above  into  equation  and  simplify  to  obtain: 

d  fs  (Ip  K4  £2  +  I,  K3  £  +  I7)1  -  3  C  Ii  (90) 

d5  {  (K3  +  K*  X?  j  (Ks  +  K„  5) 


Now  V  is  a  function  of  ' x'  in  this  case  and  is  given  by: 
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V  =  V0  cos(x/R)  +  U  sin(x/R) 


(91  ) 


Since  x  is  much  smaller  than  R,  this  may  be  approximated  as: 


V  =  V0  +  U  (x/R) 


(92) 


Substitute  into  equation  (90): 


d  fi  (Is  k4  c2  +  Is  K3  g  +  I7)l  Is  g  ((Vo/ U)  +  (e/R)  g)  (93) 

dC  {  (K3  +  K4  W  ]"  (K3  +  K„  5) 


Perform  the  integration 
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Is  (K3  +  k4  a 


K4(Ig  K4  £*  +  Ie  K3  C  +  I?) 


Vo 

U 


-3ln~ 

k4 


J 


+  ^-2  (E  £n5  -5+1) 


+  C 


(94) 


Where : 


=  ■  [fn 


+  1 


(95) 


Evaluate  the  constant: 


at  5=0,  6  =  0  =>  C  =  0 


(96) 


ft  ■  f 
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Final  1 y : 


5  = 


U  (K3  +  K„  c)2 

(Is  K„  52+  Ie  K3  l  +  I7) 


Vo 

U 


V 


(97) 


M  —  .  KsV 


-3£nE  +  -f-z( 

l\4  J\4 


£nE  -  E  +  1) 


To  evaluate  the  constants  take  as  a  first  approximation  the 
s i mi  1 ar i ty  prof i 1 e : 


fi  =  9i 


-0.693n2 

e 


(98) 


Thus 


IG  =  2.266 


(99) 


I7  =  1.133 


MOO) 


Now  from  the  experimental  results 


K3  =  0.0822 


MOD 
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K4  =  0.00439 


(  102) 


2 . 2  Comp  a  r  i  son  of  Theory  to  E  xper  i  merit 

There  is  much  experimental  data  concerning  the 
deflection  of  jets  by  crossflows.  However,  there  is  little 
data  on  the  particular  region  or  type  of  jet  that  was 
analysed  above.  At  present  the  data  consists  of  a  very  few 
runs  at  high  velocity  ratios  and  perhaps  one  or  two  initial 
points  from  low  velocity  ratio  runs.  Data  on  compound  jets 
and  wall  jets  is  even  scarcer. 

Figure  (40)  shows  a  comparison  of  equation  (78)  to 
compiled  data  of  Pratte  and  Baines  (1967).  The  agreement  is 
good  up  to  a  value  of 

V  x_  .  (103) 

U  b0 


Figure  (41)  compares  the  same  equation  to  one  profile  of 
Rajaratnam  and  Gangadrahiah  and  indicates  good  agreement  up 
to  the  same  limit  as  above. 


Finally,  equation  (97)  may  be  compared  to  the 
experimental  data  obtained  previosly  in  this  study.  This 
plot  is  displayed  on  figure  (42).  The  magnitude  of  V  was 
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Figure  42  -  Prediction  of  Jet  Deflection  for  Various  R/0  and 
Comparison  With  Experimental  Observations 
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calculated  from  equation  (3).  The  theory  seems  to  indicate 
the  magnitude  of  the  deflection  fairly  well  but  the  large 
scatter  in  the  data  prevents  any  definate  conclusions.  The 


scatter  is  believed  to  be  caused  by  difficulty  in  aligning 


the  jet  nozzle  with  the  coordinate  system. 

2 . 3  D i scuss i on  of  Theory 

The  use  of  a  moment  of  momentum  equation  to  calculate 
the  deflection  is  natural  since  the  method  is  analagous  to 
the  calculation  of  the  mean  of  a  statistical  distribution. 
The  location  of  the  maximum  velocity  point  is,  in  fact,  the 
'mean'  location  of  the  jet. 

The  physical  i nterpretaion  of  these  results  may  be 
easily  seen  by  simplification  of  the  final  differential 
equation.  Equation  (26)  for  example  may  be  simplified  by 
noting  that  the  corresponding  momentum  equation  reduces  to: 


(  104  ) 


Thus : 


dd  =  JLl  V 
dx  1 2  u 


m 


(105) 


The  deflection  is  seen  to  be  a  consequence  of  simple 


' 
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convection.  The  constant  factor  accounts  for  the  non-uniform 
velocity  distribution  and  thus  equation  (105)  states  that 
the  deflection  is  a  consequence  of  some  type  of  average 
convection.  This  equation  may  be  contrasted  with  equation 
(2).  The  essential  difference  is  that  the  former  will  apply 
in  cases  where  the  jet  characteristics  dominate  and  the 
latter  to  cases  where  they  do  not. 


« 
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CONCLUSION 

The  present  work  on  the  deflection  of  jets  by  weak 
crossflows  involved  an  experimental  and  analytical  study. 

The  experimental  results  provided  the  qualitative 
information  about  similarity  and  scale  variation  which 
allowed  the  simplification  of  the  analysis.  The  experiments 
did  not,  however,  provide  a  validation  for  the  theory  due  to 
a  large  scatter  in  the  deflections  observed.  This  scatter  is 
believed  to  have  been  caused  by  a  crude  jet  alignment 
procedure . 

The  theoretical  analysis,  involving  the  use  of  an 
integrated  moment  of  momentum  equation,  provides  a  simple, 
direct,  method  of  calculating  the  deflection  of  a  jet.  The 
method  is  only  valid  in  the  initial  region  of  a  weakly 
deflected  jet  but  there  are  numerous  practical  situations 
where  this  region  is  of  interest.  Advantages  of  this  method 
are  that:  it  is  direct,  allowing  the  use  of  solutions  for 
undeflected  cases:  it  does  not  require  an  empirical  constant 
or  function:  and  it  is  easily  applied  in  a  wide  variety  of 
cases . 


As  suggestions  for  further  study  the  most  important 
would  be  accurate  and  purposeful  experimental  work  on  the 
weakly  deflected  jet  both  to  provide  justification,  and  the 
limits  thereof,  for  the  theoretical  simplifications  for  a 
wider  variety  of  situations  and  to  validate  the  theoretical 
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calculations.  Also  of  interest  would  be  to  attempt  to  apply 
the  theory  to  problems  involving  bouyancy,  both  deflection 
of  a  bouyant  plume  by  a  crossflow  and  the  deflection  of  a 
jet  by  a  density  gradient. 
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Appendix:  Equations  of  Motion 

The  equations  of  motion  for  a  turbulent  fluid  are  (using 
tensor  notation  for  brevity): 


3u . 

1 

9t 


+  u. 
J 


I  iP_  +  I  J  + 
p  9x.  p  9x-  9x- 

'  U  1 


( A  1  ) 


9u 


i 


=  0 


(  A2  ) 


where  is  composed  of  a  viscous  and  a  turbulent  cmponent : 


=  -p  u'.u'.  +  V 
"1  J 


(A3) 


It  may  be  shown  by  order  of  magnitude  analysis  that  for 
steady  and  slender  or  boundary  layer  type  flows,  the  above 
equations  reduce  to(Pani  (1972)) 
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U2  =  .1  IE  +  I 
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(A6) 


3u  3v  .  8w  n 
3x  3y  3z  U 


(A7) 


where  the  equations  have  been  written  out  in  a  locally 
cylindrical  coordinate  system. 

This  anaysis  applies  to  crossflow  situations  where  V  is 
an  order  of  magnitude  smaller  than  the  longitudinal  velocity 
scale  and  r  i s  of  the  same  order  of  magnitude  as  the 
longitudinal  length  scale.  Equation  ( A6 )  is  simply  the 
hydrostatic  equation  and  may  be  discarded  for  the  present 
purposes . 

Integrate  equation  (A5)  with  respect  to  to  y  from  y  to 
infinity  to  obtain: 


oo 


( A8) 


Differentiate  with  respect  to  x: 
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( A9 ) 
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Thus 
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Substitute  into  equation  (A4) 
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Applying  order  of  magnitude  analysis  again,  it  is  found  that 
the  fifth  and  sixth  terms  are  negligable  and  may  be  omitted. 
Further  assume  that  the  imposed  pressure  gradient  is 
neg 1 i gable . 
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3y 


3z 


(  A12) 


3ii  +  3y+3w  (  A 1  3  ) 

3x  3y  3z  u 


For  the  the  present  purpose  it  is  useful  to  separate  out  the 
ambient  flow  in  these  equations  so  define: 


v1  =  v  -  V 


( A14) 
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u'  =  u  -  U 


(  A 1  5  ) 


where  U  is  the  ambient  co- flowing  velocity  and  V  is  the 
ambient  crossflow  velocity.  Noting  that  V  is  not  a  function 
of  y  and  U  is  not  a  function  of  x  equations  (A13)  and  ( A 1 4 ) 
may  be  substituted  into  (All)  and  (A12)  to  give: 

(U'+  U)  M'  +  (V'+  V)  M'  +  w^'  =  (A16) 

^  '  3x  v  ;  3y  3z  3y  3z 


£u'+  3v'+  3w 1 _  o  (  A 1 7  ) 

3x  3y  3z 


These  equations  and  simplifications  thereof  form  the  basis 
of  the  analytical  work  in  this  study.  A  further  conceptual 
approximation  is  made  in  that  the  these  equations  are 
applied  in  a  Cartesian  coordinate  system  whereas  they  were 
derived  in  a  locally  cylindrical  system.  Since  there  are  no 
curvature  effects  the  notation  is  identical. 


' 


